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ABSTRACT 


The dynamics of two-dimensional interacting "line" vortices is 
identical to that of the two-dimensional electrostatic guiding center 
plasma- Both are Hamiltonian systems and are therefore susceptible 
to statistical mechanical treatments. The predictions of the 
microcanonical ensemble are explored for this system. Interest 
focuses primarily on the regime of total positive interaction energy, 
which should be above the Onsager "negative temperature" threshold. 
Calculations of the probability distribution for a component by means 
of the central limit theorem are carried out in the manner of Khinchin. 
The probability distribution of a component reduces to the usual Gibbs 
distribution in the regime of positive temperatures, and is still 
explicitly calculable for negative temperatures. The negative temper- 
ature states are neither quiescent nor spatially uniform. Expressions 
for the temperature are explicitly provided in terms of the total 
particle energy and particle number. An extensive "thermodynamic 
limit" appears not to exist for the system. A BBGKY hierarchy can 
be derived for both temperature regimes. A Vlasov approximation on 

the first pair of BBGKY equations leads to the differential equation 

2 

(two dimensions) V Y + sinh Y = 0. The charge (vortex) density has 

2 

an expectation value proportional to v Y. Numerical simulations 
involving solutions of the equations of motion of 4008 particles are 


presented. 
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I. INTRODUCTION 


In 1949, Onsager made what was apparently the first application 
of statistical mechanics to a number of interacting discrete "line" 
vortices. 1 He observed that the finiteness of the total phase space 
available had peculiar implications, if one attempted to apply the 
standard equilibrium formulas. In particular, above certain values 
of the total interaction energy (a constant of the motion), the 
thermodynamic temperature was observed to be formally negative ■ 3 

Discrete line vortices as a model for two-dimensional hydro- 
dynamic flows then received relatively little attention for the follow- 
ing twenty years. Two-dimensional continuum models of hydrodynamic 
turbulence, in particular the two-dimensional Navier -Stokes equation, 
occupied the attention of many investigators, 3-10 who assumed (perhaps 
correctly) a greater degree of realism for continuum models. Even now, 
the relation between discrete and continuous models of two-dimensional 
hydrodynamic fluids remains unclear. But the subject matter of this 
article is primarily the discrete vortex model of Onsager 1 (see also 
Lin 11 ). 

The dynamics of the Onsager model has recently been shown to be 
identical to that of a two-dimensional guiding -center plasma. In this 
model, 13 very long electrostatically interacting rods remain aligned 
parallel to an external magnetic field aqd move across it with the 



4 


self-consistent if x ^ drift. This model has been the point of origin 
for a number of analytical and/or numerical calculations. 12 20 It is 
not our purpose to give a general review or assessment of these here. 

We shall be concerned primarily with a semi -rigorous applica- 
tion of the methods of equilibrium statistical mechanics to the 
guiding center plasma (or equivalently, the line vortex system). An 
earlier paper 21 presented numerical simulations of the motions of 
the system and obtained seme of the present results heuristically. 

Specifically, we are concerned with adapting the apparatus of 
the theory of probability in the form given by Khinchin 32 to obtain 
ensemble -average predictions for the states of the guiding center 
plasma at arbitrary values of the ratio of the total interaction 
energy to total particle number. We anticipate that above certain 
values of this ratio, 1 macroscopically stationary and spatially 
uniform thermodynamic equilibria do not exist . 2 > 17 The novel point is 
that, nevertheless, predictions can be extracted concerning the nature 
of the macroscopic states. It should be remarked that an understanding 
of the material in Chaps. 2, 4 , and 5 of Khinchin 22 is mandatory for a 
satisfactory understanding of the present paper. 

Section II details the dynamics of the model. Section III 
formulates the problem in the language of the theory of probability. 
Section IV applies the central limit theorem to calculate the analogue 
of the us ual Gibbs distribution for the guiding center plasma, though 
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the content of the result is radically different. Section V derives an 
analogue of the BBGKY hierarchy which applies to negative temperature 
states. Section V also treats the hierarchy in the "Vlasov limit," 
leading to a nonlinear partial differential equation which we believe 
to be of central importance in the subject. Section VI illustrates 
the results with a numerical simulation of the system. 



II. EQUATIONS OF MOTION 


The equations of motion of the line vortex system are 

die. 

1 

dt 

where x^ = (x^,y^) is the two-component vector which locates the ith 
vortex of strength I?. = K. e . x . * = x . - x . . Eqs. (l) also 

1 1 Z lj 1 J 

represent the equations of motion of the two-dimensional guiding- 
center plasma, if we make the replacement 


-E 


t. 

JL 

2 TT 





(i) 



hvc e .1? 

*L_ 

p 


( 2 ) 


where ^ is a uniform magnetic field ^ = B e , c is the speed of light, 

e. is the charge on the jth rod, and the rods, all aligned parallel 
J 

to the z-axis, have a very long length l. We are thus free to discuss 
the two systems interchangeably. In most of our discussion, we shall 
use the plasma vocabulary. 

Eq. (l) can be given a Hamiltonian form by defining canonical 


coordinates 
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and Hamiltonian function 



K.K. ton 
i J 



( 4 ) 


Hamilton' s equations 


_ 8H. a ^i . . JH 
dt ~ dp.^ ’ dt Sq^ 


( 5 ) 


then reproduce Eqs. (l). Since H is not explicitly a function of 
time, it is a constant of the motion for an isolated system. It is 
convenient to work with the total Coulomb interaction energy, which 
is proportional to H: 

i<3 i<j 


x. . 
id 


(6) 


The sum runs over all pairs of charges. Hereafter, we shall assume 

2 

there are N of +e and N more of -e in a large volume V = L . 

It will be convenient for later purposes to expand the electric 

field due to the jth charge (say) in a Fourier series, assuming 

2 

periodic boundary conditions over a box of large volume V = L : 



and the total electric field due to all charges is 


lf(x) = £. lf(x;x.) . (8) 

D J 

The integral l / E? 2 (x)dx/8rr can be interpreted as the total 

J V 

energy t if we remember to subtract off the Coulomb self-energy of 
each charge. This is the total number of charges times 

* /■ (***,» 2 s ~ 2 j w - Xt 1/k " ' (9) 

J v k 

The total Coulomb self -energy of all 2N charges (W of the +e and N 
more of -e) is, therefore, 

= UnNe 2 y , /,2 
e self _ TV L-> ' 

k 

It may be useful for some purposes to render 8 gelf finite by cutting 

off the at fiM = k , say, where k is some maximum wave 
k » ■ max. max. 

number the reciprocal of which is much less than the mean interparticle 
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2 — -* 

spacing [L /2N] s . The minimum value of k in Eq. (9) is, of course, 
2tt/L. The complete expression for the interaction energy £ is, then, 
alternatively written as Eq. (6) or 



where the E-* is again understood to be cut off at k when necessary. 
k max « 

This should be understood as a modification of the Coulomb potential 
at distances much less than an interparticle spacing. It is to be 
hoped that the essential physics of the situation is unaffected by 
this modification, and we shall not inquire into the question of 
possible collapse of the system at large negative values of the inter- 
action energy, 23 ~ 26 which might result from treating the interactions 
of positives and negatives without a cutoff. 
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III. FORMULATION OF THE STATISTICAL PROBLEM 


The equation £ = const, defines a surface of constant energy- 
in the phase space of the system. A point on this surface is con- 
veniently labelled by giving the values of the coordinates (x^,x,-,, • • • ) 
for all the charges (vortices) of the system. 

2 

If the system occupies a large but finite spatial volume V = L , 
the total volume of the phase space available to 2N charges varies as 
V . It is the finiteness of this total available phase space which 
gives the system its novel features. (By contrast, the total phase 
space available to an ideal gas in a finite spatial volume is in- 


finite. 22 ) 


If we define f ( g. ) = I die. dx, 

-W ‘ 


as 


of phase space with energy < g , then the structure 


the total volume 
function Q( P, q ), 


n( e 0 ) 



- e) 


will necessarily pass through a maximum (at g = g^, say) as g Q 
increases. The familiar classical formulas for the entropy and 
temperature. 
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S = k 0n P, Q ) + const., 

(kT )" 1 = n'(e o )/n(R Q ) , 

will predict formally that when p > the temperature T will be 
negative . 1 * 2 (h is Boltzmann's constant.) 

In this article, we are mainly concerned with the following 
question. If one makes the hypothesis of equal a priori probabilities, 
and assumes a uniform probability distribution for the phase space 
variables x^,x_, . . . over the energy shell at p - p. , what is the 
probability distribution law for a small component of the system in 
its own space? 

This is the question the answer to which in the classical cases 22 

leads to the Gibbs distribution: the well-known derivation of the 

canonical ensemble from the microcanonical ensemble. For £ < g , we 

o for 

may also expect a Gibbs distribution to result. The more interesting 
question is, What will result when P, > g^? 

Components are defined technically as groups of charges which 
interact only weakly with charges in other groups . 23 That is, the 
total energy p can be written to a good approximation as 

e = &l + e 2 + . . . + e n (12) 

in a case where there are n components, is the sum of the pairwise 
potential energies for all the charges in the kth component. The 
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microcanonical ensemble can then be defined for a system which is 
idealized as consisting of n large spatial volumes containing many 
charges each, but spatially well separated so that the charges in each 
volume interact strongly only among themselves. 

The random vector x^,Xg, ... is then assigned a probability 
distribution 


PCx^Xg, . . .) = 


6(e 0 - e) 


J^l ^2 




5(e - e) 


(13) 


It is shown by Khinchin, 23 and the proof can be readily adapted for 
this system, that the combination law for the structure functions of n 
components is 


a(s o ) - ff TT VVnW'o - £ 0 • (l4) 

V k=i J \ k=i J 

Here is the structure function for component k of energy p 

Any sum of components satisfying Eq. (12) is also a component. 
If the total system is divided into just two components, the proba- 
bility distribution law for component 1 (say) in its own phase space 



p(i) = n 2 (e o - e 1 )Me 0 ) • 


(15) 


The main computational problem we have is in evaluating asymptotic 
forms for the terms in Eq. (15), and this is addressed in the next 


section. 
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IV. APPLICATION OF THE CENTBAL LIMIT THEOREM 


The structure function fi( p) for our system is normalizable, 

in the sense that its moments can be assumed to exist. Random-phase 

estimates of the first few of these are derived in the Appendix. 

p ) is non-negative and obeys the combination law (l4), and so 

may be treated formally as a probability distribution law for a sum 

of uncorrelated random quantities P^, P,^, . . . g. , upon dividing 
2N 

each a by V . 2N is the total number of charges in the kth 

it it 

component. 2N = 2N k is the total number of charges in all com- 
ponents. 

It is important to realize that not much physical significance 

2N k 

attaches to the probability distribution laws o/V and 
£ 2N 

fi( g^)/v - ft(g Q )/v N directly. The physically significant dis- 

tribution laws are Eqs. (13) and (15). Treating 


n(e. Q ) ^(F a + P 2 + • • • + P n ) 






as the probability distribution law 


of a large number n of uncorrelated random variables is mainly a 
mathematical convenience. This convenience lies in the applicability 
of the central limit theorem from the theory of probability , 23 which 
enables one to write down asymptotic distribution laws for such sums 
as the number n of terms in the sum becomes large compared to unity. 
These asymptotic forms can then be used to evaluate Eq. (15). 
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Since our ft^'s are normalizable, our task is in some respects 

simpler than Khinchin' s. It is the lack of normalizability of fi^'s 

for systems having a classical unbounded momentum space which 

necessitates 22 introduction of the "generating functions" and 

normalizable "conjugate distribution laws" to which the central 

limit theorem can be applied. We need not make this detour, but can 

2N fc 

apply the central limit theorem directly to the 

The essential content of the central limit theorem 22 is that 
for relatively mild restrictions on the fi, 

a(e ) exp(- (e - A ) 2 /2B 1 

► T (16) 

V 21 * n-»oo (2 tt B r ) 2 

up to terms of 0(n -1 ). 2N is the total number of particles in all 
the components: 2N = 2N^ . is defined by 


k=l k=l 




2N, 


V 


and is the sum of the expectations of the first n random quantities. 

B is defined by 
n 


n n 

B n=Z b * = Z J (e K-V 

k=l k=l 


2 


2N, 


V 



1 6 


and is the sum of their dispersions. This result is independent of 
the detailed nature of the ft, . 

In the Appendix, expressions for a^. and b^ are derived within 
the framework of the random phase approximation. It is shown there 
that a^ = 0, and (see Eq. (A-10)) b^ = (N^e 2 / £) 2 /tt , where 2N^ is 
the total number of charges in the kth component. The expression 
(l6) then implies that up to terms of 0(l/n), 


"(ftp) _ exp(- 
v2H 


(17) 


In the evaluation of formula (15), we treat component 1 as a finite 
system, with number of degrees of freedom and total energy as of 
0(l). We treat both component two (all the rest of the total system) 
and the total system as very large, with number of components and 
total energy as of 0(n), n » 1. The sums E b^ and E^_g ^ are 
both 0(n) quantities which differ only by terms of 0(l), so their 
difference can be neglected. We find from Eq. (15) that 


f * \ 

B = L \ - 0(n) 

■too. / 


p(l) = 


■L exp{- (e o - g 1 ) 2 /2B} 


V 2W 1 exp{- e Q 2 /2B} 


(18) 
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As long as £^ « £ , Eq. (18) reduces to 

p (i) = "ST exp f e 0 e i/ B ' ! * ( x 9) 

V 1 

When £ q /B, which is proportional to the total interaction energy, 
is negative, Eq. (19) is the canonical distribution of Gibbs, with 
positive "temperature" ( kT) - £ q /B. The explicit value of B is 

B = \ (N k e 2 /4) 2 /TT , 

( components ) 

where 2N^ is the total number of charges in the kth component. 

For g Q > 0, the canonical ensemble again formally appears, but 
the "temperature", -B /£ , is negative. P(l) is still normalizable, 
however, because of the decay of at large £^ , for all values 

of £ q . It is meaningful to inquire into the values of ensemble 
averages computed with the probability distributions (18) or (19) 
just as it is for the more conventional Gibbs distribution. It should 
come as no surprise that some of the ensemble averaged quantities have 
different qualitative features than they do in the positive -temperature 
case, however. 

We note in passing that Taylor's 18 random phase result, £ = 0, 
(corrected for the Coulomb self-energy) has been confirmed by Eq. (19). 
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The randan phase approximation has been used in the Appendix to 

provide a specific numerical value for B, and to show that the a^ = 0. 

The canonical distribution law of a component, Eq. (19), with 

£ q /B > 0, apparently does not imply an extensive thermodynamic limit 

for the component in the conventional sense . 20 For the expectation 

value of the energy of component 1, given by using Eq. (19) in 

Jd& L p(i)a 1 (g 1 )e 1 / J d&jP(i )fl 1 ( g 1 ), can be shown via the random phase 

2 

approximation to vary proportionally to rather than N^, at least 
for small, fixed £ q /B. ( 2N^ is the number of charges in component 1. ) 

This emphasizes the necessity of rethinking from the beginning any 
macroscopic thermodynamic concepts one may wish to apply to this 
model. Variables such as the total energy are not extensive in the 
usual thermodynamic sense. 
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V. BBGKY HIERARCHY 


The canonical distribution (19) can be used as the basis for a 
BBGKY hierarchy even though £ q /B may be positive. Define the pro- 
bability distribution of the first positive charges and s e negative 
charges by, as usual, 13,20 



= V 


s ,+s 
l e 


[.•* i 
dx L , . 

J s.+l 



. -»e 

dx . . 
s +1 
e 


die® P(l) . 
1N 1 


( 20 ) 

Here, x^ is the vector position of the j0th"ion", and x^ that of the 
l th "electron" . P(l) is given by Eq. (19)* Setting 9 = - b/ , 


bn, 


s i ,s e \’ a e ^V s e } 


^d 


e 

^d 


n ° Th ■* i 
= “ T J ^a.+l 

^ i 


, ei/-»e -»i \ 

bep (x - X ) 


J s i +! J 


n 


v 

**d 


s.+l,s 
l * e 


- f / 


s +1 
e 


* ee/-»e -*e % 

bt f (x j - 
bx. e 


n 


s . , s +1 
l* e 


( 21 ) 
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Here j is any number between 1 and s g , and n Q = N^/v. 

ii 00 0i 3_0 

cp =C p = - c p = - c p is the two-body Coulomb interaction of 

Eq. (6). §(s^,s e ) is the sum of the pairwise Coulomb interactions 

among the first s^ ions and the first s g electrons. A similar hierarchy 

can be obtained by taking gradients with respect to ion coordinates. 

The hierarchy (2l) is identical with the conventional hierarchy 20 except 

that 9 may be negative, and is so in the case of most interest. 

For negative temperatures, Eqs. (2l) will not be tractable by 

the usual methods. We do not expect n^ q and n Q ^ to be spatially 

uniform for 0 < 0. It is interesting to try a "Vlasov approximation" 

on (21), approximating n^ = n. a. , n . n , n . , 

’ 2, o l,o 1 , 0 * o,2 o,l o,l ’ 

and n, . = n. n .. . The result is 

J.,1 1,0 o,l 


bn 




n 


bx-. 


T 1 ^ 


bcp 


ei/ -»e 


( -»e 
(x i " 


X 1 




^,o (x l )n o,l (x l } 



x ee/ -»e 
bcp {x 1 - 


fix-, 




bn l,o 


x -*l 

bx 1 


n o r.^e 
TJ to i 


. ei , -4i -4e s 

bcp (x 1 - x L ) 

b *-, 1 


)n o,l (i? l ) 


n ° r,-i 

" T J dX 2 


- 11 / -»1 -+1 s 

b<p(Xi-x ) 

a ”1,0^ x i )n i , 0 (x 2 ) 


bx n 


( 22 ) 
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Eqs. (22) can be rewritten as 



-»e 


n , (x, 
o,l' 1 


n 

o 

T 


/' 


dx 7 


ee/ -»e -►/ 

cp (x, - x 


)[n 0 ( i ( x') 




Adding Eqs. (23) gives 


^ n o 1^ + ^ "q 0 (5>] = 0 > (24) 

5x ’ ’ 



( 25 ) 
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and c 2 are two constants of integration satisfying = n Q q ^ 0 

which are to be determined by requiring that 



(26) 


and cp(x - x') = - (2e 2 /£)0rc|x - x ' | 

Defining p(x) = n. (x) - n n (x) , Eqs. (25) can be combined 

-L. y O O y -L 

to give 



On physical grounds, we expect antisymmetry of the charge distribution 
about the center of the box, which permits us to restrict attention 
to the case c^ = c 2 = c, where 



If we then introduce a scalar potential. 



3 ?')p(x , )cb?' 


Eq. (27) becomes 



23 


2 4T,en o c r 
V $ = I exp 


(t) ' 6X5 (■ 


e_£ 

e 


(29) 


.1 ~2 22 / 

In a more compact notation, let Y = e $/ 9 and v = v /c , with 


V* - l 


/4rrn o e^ , and Eq. (29) becomes 


V 2 Y + sinh Y = 0 


(30) 


The nonlinear differential equation (30) determines the charge- 
density distributions for the negative -temperature states. It was 
derived by a quite different method in Reference 21. Easily soluble 
in one dimension, it is apparently intractable in two, and will pro- 
bably require numerical investigation. It is to be noted that if 
the sign of 9 were reversed, the sign of the sinh Y term would be 
reversed in Eq. (30), leaving us with the Debye-Poisson equation 
in two dimensions. One effect of the negative temperature is to con- 
vert its exponentially damped linearized solutions into spatially 
oscillatory ones. 81 This set of spatially oscillatory linear solu- 
tions can be fit into the box with either perfectly -conducting or 
periodic boundary conditions to give density distributions qualitatively 
similar to those in Ref. 21, but no quantitative agreement has yet been 


achieved . 



VI. SIMULATION RESULTS 


We do not yet know analytical ensemble average predictions 
for the charge densities that would follow from the solutions of 
Eq. (30), nor analytical estimates of the fluctuations to be expected 
about the ensemble average. It is perhaps instructive to present as 
an illustration the results of numerical simulations of the high-energy 
states of the system. These are expected to be a pair of macroscopic 
counter-rotating vortices. 2 17 

More details of the program by which the equations of motion 
of 4008 particles obeying the dynamics of Sec. II are presented 
elsewhere. 31 ’ 37 The charges are confined by a perfectly conducting 
boundary (for the vortex problem, the boundary is therefore a 
streamline) inside a rectangle of width one half the length. Loading 
schemes have been developed which permit arbitrary initial values to 
be assigned to the total energy. Most of the phase space available 
to the system which is consistent with very high energies corresponds 
to configurations in which the positive charges crowd together in 
one half of the box while the negatives crowd together in the opposite 
half. The macroscopic flow pattern this corresponds to consists of 
two spatially-separated counter-rotating vortices. 

Shown in Fig. 1 is a plot of the fluid velocity vector at 
t = 0 and t = 160 (see Reference 21 for a discussion of the units) 
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for the case in which the total interaction energy for the particle - 
in-cell model is 2.5 times the self -energy. In Fig. 2, similar plots 
are given for t = 0 and t = 375 for an interaction energy 0.43 times 
the self-energy. The degree to which the macroscopic counter- 
rotating vortices form and the speed with which they form both 
qualitatively increase as the interaction energy increases. The 
capacity of such plots as Figs. 1 and 2 to display any such macro- 
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APPENDIX: RANDOM PHASE APPROXIMATION TO 0(fi Q ) 

In this Appendix, we elaborate and extend Taylor's approximate 18 
calculation of 

■0U o ) = K dx 2 - . . 6 (£ q - s) ( A-l) 

where £ is given by Eq. (10). In addition, we are interested in 

determining the moments of o(£ q ) and the character of its approach 

to zero at large £ . 

o 

The essence of the random phase approximation as applied to 

this problem is the introduction of the variables r , 

k 


p = E exp( - i is • x ) = r exp(icp ) 
k 0 W J k k 


(A-2) 


as independent variables in place of the particle coordinates x . . 

tJ 

It is assumed that convergence can be achieved by taking the number 
independent i?'s as equal to the number of degrees of freedom (in 
this case 4 n) and taking both very large. This large number can be 
identified with the number of ' s in the sums of Sec. II. The random 
phase approximation to the energy amounts to a neglect of terms with 
1c + it' ^ 0 in the double sum 


cs 



This reduces if 2 


to the random phase expression 


if 2 = 16 tt 2 r^ 2 /k 2 

k k 

independently of the values of the phases cp • 

£ 

It is then possible to show that the transformation 



and any integral of the form 


J 


dx 1 dx 2' 


(A-4) 


(A-5) 


. can be converted to one 
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/ 2 2 

J dr dr . . . , provided that the integrand depends only 
v x? 


2 

on the r ' s and not the m ' s . 


With these remarks, and in view of (A4) and Eq. (10), we see 
that we can write the expression (Al) as 


n(& Q ) = j | - f dx L dx 2 . . . exp i\(ft Q - g) 

*” GO J 


■c 


p 9 i\2n jfcVr 

2 n aht i 

r V 2 exp 2 

k 2Ke^ k 


iXUnNe 2 /v 2 2 


2 2 r -» / eXp iX£ o 

m. 2Ne k f 


All the integrals over the dr can be done, and Eq. (A -6) can be 

£ 

written compactly as 


a(e o ) 


= ^ LI dX eXP ( U ( £ ° + ^ 5 ?) + AU) ] ’ 


where 


A(X) s - Z Qm. (1 + ia 2 \/k 2 ) , a 2 s 4TiNe 2 / j0V 
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A steepest -descent approximation to Eq. (A-7) suggests itself, 

since the coefficient of X is, in the case of interest, a very large 

quantity when N » 1. This can be used to show that dfif t )/d ?, = 0 

o' o 

when = 0. Thus the threshold for negative temperatures is p = 0. 
(This is Taylor's result, 18 when the Coulomb self -energy is taken 
properly into account. ) 

Moments of the type 

< £ o"> 3 / de o e o mfi < S o> 
can also be evaluated by noting that 



(A— 8) 

In particular, (j °) = V 2 ^, and 

<e 0 l > = 0 


9 


(A-9) 
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( A-10) 


(These are the results needed in Sec. IV.) 

By writing Eq. (A-7) in the form (a 2 = 4rrNe 2 /j£V) 



(A-ll) 


We can determine the way in which £ ) falls off at large 6 q . 

(e o + £ self is never negative.) As & q ► + <*>, we may push the contour 

of X integration up in the complex plane and encounter an infinite 
sequence of (multiple) poles at X - ik^/a 2 . The closest ones of these 
are at | £ j = 2rr/L, and give an asymptotic decay of ft( & Q ) as 
P(e o ) exp £ o 4TT 2 /a 2 L 2 J , as e o * oo , where p( t Q ) is a poly- 

nomial of degree at most equal to the multiplicity of the pole at 
X = i(2TT) 2 /L 2 a 2 . For a range of g, > 0 this means that Eq. (19) is 
normalizable directly, but above that range of t Q we must use (l8) 
to insure normalizability. 
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FIGURE CAPTIONS 


The fluid velocity vector field as a function of spatial 
position for the case of 4008 particles with interaction 
energy =2.3 times the self energy: (a) initially and 

(b) at t = 160.0. The velocity field at all later times 
is qualitatively similar to that at t = 160. The lengths 
of the arrows are normalized to V-MAX, the maximum dimen- 
sionless velocity present at any instant. This case should 
correspond to negative temperatures. 

The fluid velocity vector field as a function of spatial 
position for the case of 4008 particles with interaction 
energy = 0.43 times the self energy: (a) initially and 

(b) at t = 375- The lengths of the arrows are normalized 
to V-MAX, the maximum dimensionless velocity present at 
any instant. This case should correspond to positive 
temperatures. The velocity field in (2b), which is similar 
at all later times, differs from that in (lb) by the ab- 
sence of the two larger vortices. 
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